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HgTe quantum wells and surfaces of three-dimensional topological insulators support Dirac 
fermions with a single-valley band dispersion. In the presence of disorder they experience weak 
antilocalization, which has been observed in recent transport experiments. In this work we conduct 
a comparative theoretical study of the weak antilocalization in HgTe quantum wells and topological 
surface states. The difference between these two single-valley systems comes from a finite band 
gap (effective Dirac mass) in HgTe quantum wells in contrast to gapless (massless) surface states in 
topological insulators. The finite effective Dirac mass implies a broken internal symmetry, leading 
to suppression of the weak antilocalization in HgTe quantum wells at times larger than certain tm , 
inversely proportional to the Dirac mass. This corresponds to the opening of a relaxation gap in 
the Cooperon difi'usion mode which we obtain from the Bethe-Salpeter equation including relevant 
spin degrees of freedom. We demonstrate that the relaxation gap exhibits an interesting nonmono- 
tonic dependence on both carrier density and band gap, vanishing at a certain combination of these 
parameters. The weak-antilocalization conductivity refiects this nonmonotonic behavior which is 
unique to HgTe QWs and absent for topological surface states. On the other hand, the topological 
surface states exhibit specific weak-antilocalization magnetoconductivity in a parallel magnetic field 
due to their exponential decay in the bulk. 

PACS numbers: 



I. INTRODUCTION 

Recently discovered materials - graphenejii^ two- 
dimensional (2D)2i"— and three-dimensional (3D)2i"— 
topological insulators (TIs) ^'^i^^ - exhibit a Dirac-like 
band dispersion which is responsible for their unusual 
electronic and optical properties. In graphene the 
low-energy electron spectrum can be approximated by 
two spin-degenerate Dirac cones at the corners of the 
Brillouin zonei^ The 2D TIs have been realized in 
HgTe quantum wells (QWs)^i^ which have a single 
double-degenerate Dirac valley, as predicted by band- 
structure calculations and inferred from transport mea- 
surementSfiS, The double degeneracy of the HgTe QW 
bands allows for an energy gap at the Dirac point, with- 
out breaking time-reversal invariance, which paves the 
way to study Dirac fermions with a finite (positive and 
negative) effective mass and related mass disorderii2r— 
In comparison with HgTe QWs, the ideal 3D TI exhibits 
a single non-degenerate gapless Dirac cone on the surface 
of the material, whereas its bulk is insulating.^" In this 
case, the opening of the gap in the surface Dirac spectrum 
requires time-reversal symmetry breaking and has been 
predicted to cause the surface quantum Hall eSec^^— 
and rich magneto-electric phenomenap^i"— related to ax- 
ion electrodynamics.^* 

The number of the Dirac valleys is an essential fac- 
tor in understanding quantum electron transport in dis- 
ordered samples. The transport studies of graphene 
have reported weak localization^^, and more complex 
quantum-interference patterns^ instead of the antilo- 
calization effect expected for the symplectic universal- 



ity class (e.g. Refs. [3ll - l47l ). Such a situation can oc- 
cur if the two graphene's valleys are coupled by scatter- 
ing off atomically sharp defects.— i^i^i^i^^ In contrast, 
in single-valley Dirac systems such scattering processes 
are forbidden, and recent experiments on HgTe QWs^i^ 
and 3D TIe^°"^^ have found a positive (antilocalization) 
quantum- interference conductivity. 

In this work we conduct a comparative study of the 
weak antilocalization (WAL) in HgTe QWs and on sur- 
faces of 3D TIs. The goal is to elucidate the difference 
between these two systems which comes from the finite ef- 
fective Dirac mass in HgTe QWs in contrast to the mass- 
less surface states in 3D TIs. Like in the conventional 2D 
electron systems (2DESs) with Bychkov-Rashba or Dres- 
selhaus spin-orbit interactions, '^^''^^ the WAL in HgTe 
QWs and on surfaces of 3D TIs reflects the broken rota- 
tion symmetry in relevant spin space. However, in addi- 
tion to the lack of this continuous symmetry, the effec- 
tive Dirac mass in HgTe QWs implies a broken discrete 
symmetry which for a single-cone system would play the 
role of time reversal. In this sense, there is a formal anal- 
ogy between the effective Dirac mass and an out-of-plane 
Zeeman field in a 2DES4^'^"^ Therefore, by analogy with 
weak ferromagnets^^i^ we find that the WAL in HgTe 
QWs is suppressed at times larger than certain tm i in- 
versely proportional to the effective Dirac mass. Such 
suppression is however absent for the massless surface 
states in 3D TIs, which can be used to experimentally 
differentiate between the two systems. 

Before going to the calculation details given in Sec. IIIII 
and IIVI in the next section we briefly announce some of 
our results for HgTe QWs and TIs. 



2 






M= +5 meV 




M= meV 




M=-7.5 meV 




Al=-15 meV 



n(xlO'^cm-^) 





n=lxlO"cm-^ 




n=1.5xlO"cm-^ 




n=2xl0"cm-^ 




n=2.5xlO"cm-^ 



10 



Al(meV) 




FIG. 1: (Color online) Relaxation gap [in units of inverse 
life-time , see Eq. ((2)1] versus carrier density n (a) and band 
gap M (h); A = 380 meV-nm and B = 850 meV-nm^ (from 
Ref.Hl). 



II. OVERVIEW OF THE RESULTS 
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FIG. 3: (Color online) Quantum-interference conductivity 
correction Saxx [in units of e^/nh, see Eq. ((!])] versus nor- 
malized dephasing rate tq/tip for different carrier densities 
and = — 10 meV (a) and for different band gaps and 
n = 1.9 X 10" cm-^Tb); A = 380 meV-nm and B = 850 
meV-nm^ (from Ref . fig) . See also Eq. ((22)) for r and to and 
Eq. ((53) for a. 



To calculate the quantum-interference (Cooperon) con- 
ductivity correction, Saxx, we adopt the effective Dirac 
model of HgTe QWs^ and obtain the following expression 
for Saxx- 




FIG. 2: (Color online) Quantum-interference conductivity 
correction Saxx [in units of e^/ivh, see Eq. ((T])] versus car- 
rier density n and band gap M; A — 380 meV-nm, B — 850 
meV-nm^ (from Ref. (T^ and tq/t^ = 0.01 [see also Eq. ^ 
for r and tq and Eq. ((55)) for a]. 



Here the symmetry-breaking-induced relaxation gap 
is proportional to the total mass term, A4 + Bk^ , in the 
effective Dirac model of HgTe QWs.'' Ai is the band gap 
at the Dirac point, Bk^ is the positive quadratic correc- 
tion accounting for the curvature of the filled conduction 
band {k^ is the Fermi momentum determined by the 2D 
carrier concentration n), whereas constant A determines 
the linear (Dirac) part of the spectrum (r the elastic life- 
time). In Eq. ([T|) the factor of 2 accounts for the double 
degeneracy of the Dirac valley in HgTe QWs, constant 
a approaches —1/2 for M + Bk'^ — )■ (as discussed in 
detail in Sec. IIII() and is the dephasing time. 

The unique feature of the HgTe QWs is that the band 
gap M can take both positive and negative values de- 
pending on the QW thickness i^ii^ Therefore, the relax- 
ation gap tJ^ ([2|) exhibits an interesting nonmonotonic 
behavior as a function of both A4 and carrier density n 
[see Fig. [Tj , vanishing when these parameters satisfy the 
condition, 

M + 2TTnB = 0. (3) 

It represents a line in parameter space (A4 , n) on 
which conductivity ([T|) reaches the maximum Saxx — 
(c^/tt/i) ln(r<p/r) [see also Fig. [2]. Such a nonmonotonic 
behavior of Saxx{n, M) can be used to experimentally 
identify the symmetry breaking and the resulting re- 
laxation gap tJ^ . In particular, the predicted carrier- 
density dependence should hold for the QWs where the 
potential impurity scattering is much stronger than scat- 
tering off random gap fluctuations. This regime can be 
identified from the carrier-density dependence of the QW 
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mobilityiii As to the dependence on the gap M, it can 
be extracted from sample-to-sample measurements. The 
band structure of MBE grown HgTe QWs is controllable 
to a great extenti^'^'i^ For the experiment we suggest here 
one should select several QWs with distinctly different 
gaps and comparable dephasing times (inferred from the 
temperature dependence of the conductivity) and other 
band structure parameters (inferred from the Hall and 
Shubnikov-de Haas measurements). Alternatively, the 
presence of the relaxation gap tJ^ can be established 
from the dependence of Sa^x on the dephasing rate 
which is directly related to the temperature dependence 
(e.g. the dephasing rate due to electron-electron interac- 
tions is linearly proportional to the temperature) .~ The 
dependence of Saxx on is shown in Figs. [3l[a) and 
(b). In these figures the upper curves correspond to the 
logarithmically divergent Saxx with tJ^ — 0. In con- 
trast, for finite tJ^ (rest of the curves) the conductivity 
correction shows only weak dependence on the dephasing 
rate. 

As to the 3D TIs, we focus on compounds 612863 and 
Bi2Te3 where the surface states can be described by the 
effective two-band Dirac Hamiltonian, accounting for the 
hexagonal warping of the bands [see, e.g., Rcf. 55^5^. 
The warping term is cubic in momentum k and enters for- 
mally as the Dirac mass term. However, since it preserves 
the time-reversal symmetry, we find that for weak warp- 
ing the quantum-interference conductivity correction has 
the same form as for the conventional 2DES with spin- 
orbit interactions^ii^i^ 



Saxx = -a— In a = --. (4) 
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FIG. 4: (Color online) Energy bands of a HgTe quantum well 
[see Eq. and text] in meV versus in-plane wave-numbers 
kx and ky in nm~^. The Fermi level lies in the conduction 
band at _B = 0. We choose A = 380 meV-nm, B = 850 
meV-nm^ V = 700 meV-nm^ (from Ref. 16), and = 100 
meV. 

III. HGTE QUANTUM WELLS 

A. Effective Hamiltonian 

We use the effective 4-band Hamiltonian of HgTe wells^ 
which can be written in the following form: 

fH + H, \ 

HHsTe-[^ H + hJ' 



The specific of the surface state shows up mainly in the 
dependence of the conductivity AaxxiB) = SaxxiB) — 
S(Jxxi(^) on magnetic field B applied parallel to the TI 
surface: 

Such dependence refiects the finite penetration length, 
A, of the surface state into the bulk, i.e. the magnetic 
flux through the effective width of the surface state [see 
Eq. ([5]) for B^^ , where £tr and are the transport mean 
free path and dephasing length, respectively]. Quantum 
transport in the in-plane magnetic fields has been studied 
theoretically for disordered metallic films^ and electron 
quantum wells.-S? The present case differs from the previ- 
ous studies in that the topological surface states have a 
different microscopic profile of the transverse wave func- 
tions. We discuss the dependence of Aaxx on the mag- 
netic field orientation in Sec. II VI in connection with recent 
experiments on Bi2Te3 (Ref. i52i ). 



H = (T{Ak + Mi^z) + {Vk^ -E^)(7°, (7) 
H^-cT{Ak + Mi^z) + {Vy-E^)a'>, (8) 
H, = l/(r)aO. (9) 

The two diagonal blocks in HugTe describe pairs of 
states related to each other by time reversal symmetry 
(Kramers partners). In the upper block the Hamilto- 
nian H has a matrix 2x2 structure with Pauli matrices 
u^,v,^ and unit matrix cr*' acting in space of two lowest- 
energy subbands of the HgTe quantum welli^ an s-like 
electron one \Ei,^) and a p-like heavy hole one |). 
For the lower block the basis states have the opposite spin 
projections: i) and |). The linear term in 

H (proportional to constant A and momentum operator 
k — — iVr) originates from the hybridization of the s- 
and p-like subbands. is the effective Dirac mass: 

Mi^ = M+Bk\ (10) 

where constant M determines the band gap at the 
gamma (k = 0) point of the Brillouin zone (see Fig. |4]). 
The positive quadratic terms Bk^ and Vk^ take into ac- 
count the details of the band curvature in HgTe quantum 
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wellsj^i^ Finally, Hi in Eq. ([6]) accounts for interaction 
with randomly distributed short-range impurities. The 
impuritity potential V{y) is characterized by the correla- 
tion function, 



av-v')^{{V{v)V{v'))) 



Ck - / C(r)e-^'''-dr 



TTiVro' 



5{y-t'), (11) 
(12) 



parametrized in terms of the characteristic scattering 
time To and the density of states (DOS) at the Fermi level, 
N , for one Dirac cone [brackets ((...)) denote averaging 
over impurity positions and (^k is the Fourier transform 
of the disorder correlation function]. 

We emphasize that the mass term M.^ ([TO|l is sym- 
metric with respect to momentum inversion k — > — k. 
Hence, Hamiltonian H does not possess the symmetry 
under transformation 



k, cr — > — k, — cr. 



(13) 



Within a given block of Eq. (jS]), i.e. in subband pseu- 
dospin space, such a transformation plays the role of 
" time reversal" . At the same time, the real time-reversal 
symmetry is ensured by the matrix form of HugTe ©. 
Physically, this means that the Kramers partners reside 
on two Dirac cones superimposed at k = point. ^® Note 
that the zero off-diagonal elements in HugTe (© imply 
conservation of the spin projections of \Ei) and \Hi) 
subbands, which is a good approximation for symmetric 
HgTe wells In this case, each of the Dirac cones 

contributes independently to transport processes, which 
is accounted for by the factor of 2 in the expressions for 
the conductivity corrections discussed in subsection llll Dl 



B. Disorder-averaged single-particle Green's 
functions and elastic life-time 
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FIG. 5: (Color online) Diagrammatic representations for (a) 
equation for disorder-averaged Green's function (thick line) 
in self-consistent Born approximation (thin line is the unper- 
turbed Green's function of the disorder-free system, dashed 
Une is the disorder correlator), (b) bare and dressed Hikami 
boxes for the Cooperon correction to Drude conductivity, (c) 
Bethe-Salpeter equation for the Cooperon, and (d) equation 
for the renormalized current vertex in the ladder approxima- 
tion. 



Here the Green's function G, 



R/A 

Ok 



describes a conduction- 



band carrier with dispersion = \/ A^\f? 
Ep in the absence of disorder [index s = ±1 labels the 
Kramers partners residing on the different Dirac cones]. 
The valence band contribution is neglected under as- 
sumption that the energy separation between the valence 
and conduction bands is much bigger than the disorder- 
induced band smearing, 



We begin by calculating the disorder-averaged re- 
tarded/advanced Green's functions G ' for an n-type 
HgTe well under weak-scattering condition 



kpVpT > 1, 



(14) 



where r is the elastic scattering time and and kp = 
^/2im are the Fermi velocity and wave-vector, respec- 
tively. In the standard self-consistent Born approxima- 
tion (see, e.g. Refs. IgiHb^ ) the equation for G^^"* is 
shown diagrammatically in Fig. [SJa). In k representa- 
tion it reads 



Gy, 



G, 



Ok 



G, 



Ok 



^R/A RIA 

,^ sn Ck-k'Gk/ Gk , (15j 



Ml 



2{Ep-Vkl):^h/T. (17) 



Because of the large band-structure constant A ~ 380 
meV-nm [see, e.g. Ref. [l6|] the requirement (flT)) can be 
satisfied simultaneously with the weak-scattering condi- 
tion ([H]) . We also note that the matrix structure of Gg^^ 
(|16p accounts for the carrier chirality and is of primary 
importance throughout the paper. 

The solution to Eq. can be sought in the form 



Gk 

and S 



2 

R/A 



. R/A R/A 

-(^^kjgy: , where 5^ 



i(a°- 

satisfies the equation 

Ck-k' 



i/(6-a-Sk ) 



R/A 



1 -I- ekCk' dk' 



Ck 



^R/A 



(18) 



G, 



R/A 

Ok 



1 cr 
2 







o-ek 



-,ek = s 



A\l + Mwz 

^A^k-^ + Mi' 



|ek| 



1,(16) 



Approximate solutio: 



,61.62 



of the latter equation near the 



Fermi surface, |k| « k^, yields the disorder-averaged 



5 



Green's function as 



2 



"R/A 



where r is the elastic hfe-time given by 

1 2^ _R/A N 



-1' 



1 + ej_„ • e^„/ 



ih 
Zt 



(20) 



II n "=11 n' 



sni/l 



(21) 



In Eq. (|20|) the unit vectors n and n' specify the direc- 
tions of the incident and scattered momentum states, re- 
spectively, and e..n and ej^n are the in- and out-of-plane 



components of the unit vector ekf 



e_Lr 



For 



isotropic CfeF |n-n'| (HSI and Mtp n (ITOl) one finds the 
elastic time 



To 



1 + 



(22) 



which is shorter than the disorder-related time scale tq 
[see, Eq. (dJ)]. This is due to the allowed backscattering 
into an opposite momentum state (nn' « — 1) which is 
absent in the gapless case.— The backscattering is the 
consequence of the symmetry breaking due to the mass 
term. The strength of the symmetry breaking is con- 
trolled by the out-of-plane component exn of the unit 
vector efc^.n [see, Eq. (HH)]. 



C. Cooperon 

The quantum-interference corrections to the Drude 
conductivity can be expressed diagrammatically by the 
Hikami boxes shown in Fig ^}p. Apart from the 
single-particle Green's functions G^J they in- 

volve the disorder-averaged two-particle Green's func- 
tion, Capa' fji (q) , known as the Cooperon. In this subsec- 
tion we will set up and solve the equation for Capa' /3' (q) • 

For the potential disorder defined by Eq. (ITT|) the dia- 
gram for the Cooperon equation (see, Fig [St) is read off 
as follows 



a^.','(q) - -Saa'Spp^ + ^ y (2^ 

^ H '^"7' + ^^)'^/35' (q - k, e)CYS'a'f3' (q), (23) 

'y' 5' 

where the Greek indices label the states in pseudospin 
(cr) space. The prefactor t^/tq in the first term is due 
to the chosen normalization of Ca^„'^'(q). To solve 
Eq. it is convenient to first expand the Cooperon 



in the orthonormal eigenfunctions of the pseudospin of a 
two-electron system: 

c^p^,p, = J2c''KpK%', E*i/3*"/^ = ^^- (24) 

ij a 13 

The indices i,j — 0,x,y,z label the pseudospin-singlet 
(0) and pseudospin-triplet {x, y, z) states. The conduc- 
tivity corrections will be eventually expressed in terms 
of the coefficients C*-' for which we derive the following 
algebraic equations (see, also, Appendix IXj) : 



C^^(q) = -*,, 

To 



(25) 



- ^ Tr((aO - rTe_fe,.„)a*(aO + ^efc,.„)a^)C^J"(q), 



4Tr 



where the square brackets stand for integral over the di- 
rections of the momentum k = fc„n on the Fermi surface: 



27r 



Q 27r 1 — itljJ -f iTVj^ n • q 



(26) 



Evaluating the traces of the products of the Pauli matri- 
ces in Eq. (|25|) we find 



To 2ro 

+ ^ E ((e+-e-+^e+xe_)•b)C^^ (27) 
2to , ^ 

b—x,y,z 

= —5a, + 7^((e+ - e_ - ie+ x e_) • a)C°^' 
2to 



To 

— y 

2Tn . ^ 



((1 + e_|- • e_)(a • b) — i(e+ + e_) • a x b 



h—x,y,z 

-(e+.a)(e_.b)-(e_.a)(e+.b))C^^-. 



(28) 



We separated the singlet C^-' and triplet C'^"''^^^ Cooper- 
ons with respect to the first index so that a, b run over 
x, y, z only. Respectively, vectors a and b run over the 
unit vector basis of the Cartesian system. We also intro- 
duced a convenient shorthand notation 



e± = e±kf 



(29) 



As discussed in subsection IIII Al the specifics of the ef- 
fective Hamiltonian for HgTe quantum wells consists in 
the symmetry of the mass term pop . Being an even func- 
tion of k, it breaks the symmetry of Hamiltonian H in 
Eq. ([7]) under transformation k, cr — > — k, — cr. The sym- 
metry breaking is encoded in the unit vectors e± (j29p 
which have antiparallel in-plane and parallel out-of-plane 
components. 



e± = ±6,1 -I- ej_. 



(30) 



where e^^ and ej^ are defined by Eq. ([?T|) . In view of the 
identities 

e+ + e_ = 2e^, e+ - e_ = 2e|^ , 

e+ • e_ = 1 — 2e^, e+ x e_ = 2e x e^, (31) 
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Eqs. ^ and ^ reduce to 

To To " 



To 



h—x,y,z 



To To'' II II 



(32) 



(33) 



(ei(a • b) - ie± • (a x b) + (ej • a)(e|| • b) 



'^0 y, 

o—x,y,z 

ie^-t,){e^-h))C'\ 



or, explicitly, 



{ey - ie^e_L)C^^ = tSoj, 



3x ZGyC^ 



= t6^ 



(34) 



(35) 



To 



(36) 



To 



(37) 



For the quantum-interference conductivity corrections we 
will only need the q- and w-dependent diagonal Cooper- 
ons (7"°, C^^ and C^^ . Each of them is obtained from 
coupled Eqs. ([M)) - where index j should be set to 
0, X or y, respectively. The coefficients in these equa- 
tions are evaluated by expanding^ the denominator in 
Eq. ()26p in the small Cooperon momentum q and fre- 
quency LU, 



Tv^n • q < 1, 



TW < 1. 



(38) 



In doing so, we keep the lowest order terms that yield the 
nonzero angle average (...) and compete with the small 
symmetry-breaking parameter ej_ [see, Eq. (1211) ]. Under 
these approximations we obtain the following expressions 
for the diagonal Cooperons: 



C"'^(q,^) 



2e2 



Dq2 



'M 



(39) 



D = Dr{2 + 5el + ei)/el 



Dr = vlT/2, (40) 



C--(q,w) 



2r 2 - e^l cos^ ( 
e2 2 + 5e2 + t 



(41) 



2t(t 



-1 

M 



2 \2 „6 „;„2 



(l + 3e^)-^-e^ sin^ 



5ei + el 



C"'^(q,^), 



C^^(q,c.) 



2t- 2 — sin^ i 
^ 2 + 5e2^ 



(1 + 3e2^)2 - cos^ 



(42) 



5ei 



■C°°(q,c.), 



where angle 0q in Eqs. (j4T|) and (|42|) indicates 
the Cooperon momentum direction: q = IqI ' 

(cos 0q, sin (/)q, 0). 

Let us analyze Eqs. p9p - (j42|) . The symmetry break- 
ing has a three-fold effect on the Cooperons. First, it 
results in a relaxation gap tJ^ in the singlet Cooperon 
C"° (|39|) . which implies suppression of the quantum in- 
terference for times larger than ryvi (even in the absence 
of the phase breaking, i.e. for cj — > 0). Second, the 
symmetry breaking affects the diffusion constant D in 
Eq. (j40|) . The diffusion constant renormalization comes 
from the off-diagonal Cooperons. In the absence of 
the symmetry breaking (i.e. for — 0) one finds^i^ 
D = 2Dr = v'^Ttr/2 with Ttr — 2t. Finally, the expres- 
sions for the triplet Cooperons C^"^ ^ and ^ 
contain additional terms oc t/tj^ = 2e^, remaining fi- 
nite in the limit w — >■ 0. Despite the smallness of the 
parameter , these terms give a noticeable contribution 
to the net conductivity correction. We will return to this 
point when discussing Eq. (|53p in the next subsection. 



D. Hikami boxes and net conductivity correction 

We now turn to the evaluation of the Hikami boxes 
for the conductivity corrections, shown in Fig. [SJd. With 
the Cooperon defined by Eq. p3|) the first and second 
diagrams in Fig.[S]3 correspond to the following analytical 
expressions: 



•M 



7 



de 



[/(e) -/(e + M] 



(27r) 



C/3^'77'(q, w) 



7r2iV2ror2 



de 
27ra; 



[/(e) -/(e + M] 



(2^ 



■C/3/3'7-y'(q,tj) 



dk 



dk' 



(27r)2 y (27r) 



where VjJ is the matrix current vertex renormahzed by 
disorder (jlll) in the usual ladder approximation (see, e.g. 
Ref. [62| and diagram in Fig. [5ji) and /(e) is the Fermi 
distribution function. We will skip the details regarding 

I 



the third diagram in Fig.[S^ since it finally gives the same 
result as Eq. (l44t . 

Evaluation of the k integrals in Eqs. and (jH]) 

yields (see, also, Appendix [B|) : 



X (ct" + cre_fc^.„0/Ji/3[('^" + o-efc^-n')Vfc^.„'(<^" +o-efc^-n')]., / 7:SlC'/3/^'77'(q)> (46) 



where the bar denotes averaging over the momentum di- 
rections n on the Fermi surface: (...) = J^^ ...d(/)n/27r. 

(2) 

We note that the correction 6axx is entirely due to the 
carrier chirality. If we omit the chirality matrix crekpn 
in Eq. (1461) , the independent averaging of the current ver- 



tices .„ gives 6ai; 



(2) 



0. 



The renormahzed vertex V^^ acquires the standard 
prefactor Tt^/T, where r^^ is the transport scattering 
time: 



1 N f 

— = ^ / d0„'(l - n • n')CfeF|n-n'| 
Ttr n J 

1 -I- ej_„ • ej_„/ + e,,a ■ e,,^. 



l + 3ei 

2to 



(47) 
(48) 



and satisfies the identity 

(t7"+crefe^.„)V^"^.„(a° + (Tefe^.„) = 2^i;^n, (49) 

X {a° + crekp.n), 
which helps to perform the averaging in Eqs. (j45l) and 
The conductivity corrections Sai]] and Sa^^^J can 



then be expressed in terms of the singlet and triplet 
Cooperons as follows 

^-^^-^-^(v)7(^[^^(^-i)^"°(-^) 

- 2n2(l - e2)C--(q) - 2^2(1 - e2)C^^(q)] ,(50) 



2e^D. 



f dq 



nh TTo "II ^'-^^ J (27r)2 
[(1 + ei)C"°(q) - C--(q) - eiC^^(q)] ,(51) 



where n2 = 1/2. We have also summed up the con- 
tributions of both Dirac cones of the HgTe QW spec- 
trum, which yields the factor of 2 in front of the integrals. 
Noticing further that on average over the directions of q 
the triplet Cooperons (|4T|) and p2)) coincide, we express 
the net conductivity correction in the form: 



^2 , 



2e^Dr /Ttr\'^ f qdq 



27r 



(52) 



1 - ei - ^ j C"0 (q) - h + ei - ^ j (q) 



8 



where (...) = J^^ ...d(j)q/2Tr with <f>q defined in text af- 
ter Eq. . Note that the terms cx 1 =F come from 
Eq. (1501) for Sa^'J , whereas the terms oc e^/2 come from 

Eq. dSI]) for (5o-E\ muhiplicd by 2. 

In Eq. ([5^ the singlet Cooperon C''°(q) results in a 
positive conductivity correction (antilocalization) com- 
ing from the pairs of states with antiparallel projections 
of (T, which prevents the constructive interference.— 
This can also be viewed as the manifestation of tt Berry 
phasei^ In contrast, the conductivity correction due to 
the triplet Cooperons is negative (localization) because 
the interference of the states with the parallel projections 
of (7 is constructive. When integrating the singlet and 
triplet contributions in Eq. (15^ we insert Eq. (15^]) and 
the terms with the diffusion pole structure (cx C°"(q)) in 
Eqs. dS]) and dHJ)^ With the upper integration cut-off 
{Dt)~^^^ and replacement —ioj in Eq. ([5^ . we 

obtain Eq. ([T|) for the logarithmic correction to the Drude 
conductivity, where the relaxation gap tJ^ is defined in 
Eq. (j39p and the prefactor a is given by 

a = ^ x(l±^V (53) 

I 2_ fj_ 1\ l + 3ei2(l + 3ei)^-e;; ] 
^[tm^tJ 2 + 5ei + el J' 

This expression is factorized into the three parts that 
have different origins: the first comes from the renormal- 
ization of the diffusion constant D in Eq. (j40|) . the second 
is due to the vertex renomalization [see, Eqs. ([52]) . (|22|) 



and (|48)l ] , and the third includes the contributions of the 
three Hikami boxes in Fig. [5}3 with the interplay of the 
singlet and triplet Cooperons [see, Eq. ([52]) ]. It should be 
noted that in the presence of the k, cr — > — k, —a symme- 
try there is partial cancellation of these three factors^, 
yielding a = —1/2 for the symplectic disorder class^ 
(see, also Eqs. (|75p and d75|) for topological insulators in 
Sec. HVl) . 

Equation ([55]) . as well as the equation for the conduc- 
tivity correction ([T}, is valid under conditions 

t/tm = 2ei < 1, t/t^ < 1, (54) 

when the carrier diffusion is limited by the time-scale, 
min{T;^, r^vf}, longer than the elastic life-time r. In 
particular, for = and t/t^ — >■ the parameter 

a — > —1/2, and we recover the result Saxx = 2x In ^ 
for the symplectic class [the factor of 2 accounts for the 
two Dirac cones, see Eq. ([6])]. For a finite the bro- 
ken k, cr — > — k, — (T symmetry leads to the deviation of 
a from —1/2 [see. Fig. [6]. The deviation is quite signifi- 
cant because the expansion in powers of involves large 
numerical coefficients: 

This behavior can be seen as the crossover between the 
symplectic and unitary classes. 



IV. SURFACE STATES IN TOPOLOGICAL 
INSULATORS 



Effective Hamiltonian 



a 



-0.5 




n(xlO'Vm-2) 



M= +5 meV 
M= meV 
M=-5 meV 
Al=-10meV 
A1=-15 meV 



a 



10 




Al(meV) 





n=2xl0"cm-^ 




n=2.5xlO"cm-^ 




n=3xlO"cm-^ 




n=3.5xlO"cm-^ 




n=4xl0"cm-^ 



FIG. 6: (Color online) Parameter a (|53|) versus carrier density 
n (a) and band gap M (b); A = 380 meV-nm, B = 850 
meV-nm^ (from Ref. 13) and tq/t^ = 0.01. 



Unlike HgTe wells the spectrum of surface states in 
topological insulators (TIs), such as Bi2Se3 or Bi2Te3, 
consists of a single Dirac cone. Consequently, the effec- 
tive Hamiltonian for the surface state in TIs^i^ has the 
form of the single-block Hamiltonian of Eq. ^ : 



Hti — H + Hi 



(56) 



where H and Hi are given by Eqs. ^ and ([9]), respec- 
tively. There are two further distinctions between Hamil- 
tonians for Bi2Se(Te)3 and HgTe wells. First, here the 
basis functions correspond to ^ and — ^ spin projections, 
Pauli (cr^'*''^) and unit (a") matrices act on real 



i.e 



spin indices. Second, the A^g.-term in Hamiltonian H 
in Eq. ([7]) is cubic (odd) in momentum k: 



55.56 



^± — kx i iky , 



(57) 



causing no gap at k = 0. This term does not break the 
k, <T — >■ — k, — cr invariance, which is now the real time- 
reversal symmetry. Instead, it causes hexagonal warping 



FIG. 8: (Color online) Polar plots of surface-state spectrum 
_iooo E{k, <l)n) (in units of energy Ak) as a function of momentum 

direction specified by angle £ (0, 2-k) [see, also, Eq. (f58)) ] . 
For W ^ the spectrum shows hexagonal warping (we chose 
W^k'^/2A^ = 0.4 and Vk/A = 1). On average over all angles 
the warping results in larger E{k,(j)-n) (dashed circle) com- 
pared to the W = Q case (solid circle). 



FIG. 7: (Color online) Energy bands of a topological surface 
state in meV [see, Eqs. (|56|) . (|57p and text] versus in-plane 
wave- numbers kx and ky in nm~^. The Fermi level lies in the 
conduction band at = 0. We choose A = 300 meV-nm, 

= 50 meV-nm^ (see, e.g., Ref. [H), 7? = meV-nm^ and 

= 600 meV. 



of the surface-state spectru m^^i^^ (see, also Figs. [7] and 



(1 + cos6(/>„) +1?^^ (58) 



We will treat the warping as weak perturbation onto the 
isotropic spectrum, assuming the smallness of the param- 
eter: 



2 1,4 



< 1. 



(59) 



Then the main effect of the warping consists in the in- 
crease of E{k,<j)n) on average over aU angles (see, 
dashed circle in Fig. [8]). This amounts to replacing 
by its angle average, 



Ml 



Ml 



2l6 



(60) 



in Eq. (|58p . In fact, the same replacement can be done in 
all even functions of A^k, e.g. Fermi momentum, Fermi 
velocity, DOS etc. Then, the specific of the surface states 
is captured by the odd carrier chirality, ere = — cret, 
in Eq. (fT6|) . Therefore, the results of the integration over 
the single-particle momenta k (given in Appendices |A] 
and|B]) apply also in the present case. This allows us to 
use Eqs. ((20|) and (|47l) for the scattering times, Cooperon 
equations ([27|) and ((28|) as well as the Hikami boxes (|45|) 
and ()46p to obtain the weak antilocalization conductivity 
corrections for the surface state in TIs. 



B. Disorder-averaged Green's functions and 
scattering times 

As in the case of HgTe QWs, we assume that the Fermi 
level lies in the conduction band of the topological sur- 
face state [Fig. [7] and the condition of weak scatter- 
ing, Eq. (ITi|) . is fulfilled. With perturbative treatment 
of the warping described above, the retarded/advanced 
Green's functions and the elastic life-time for the surface 
state in a TI are given by Eqs. (HH) - ([2T]). The time- 
reversal symmetry of the surface state is encoded in the 
odd momentum dependence of the out-of-plane vector 
e^n. There is no backscattering in this case^, and the 
elastic life-time (PU)) for isotropic Ckp n-n' is 



(61) 



In the similar way, Eq. (|47p gives the transport scattering 
time as 



Ttr = 



To 



2Tn 



l-e2/2 l+Ml^JA^kl 



(62) 



Its dependence on e^^ (and, hence, on the variance 

■^If-ii) originates from the anisotropy of the Dirac- 
fermion scattering probability in momentum space. 



C. Cooperon 

To write the Cooperon equations ([ST)) and ([SS)) for the 
surface state in a TI we note that the time-reversal sym- 
metry implies the identities: 



e = — 
e+ • e = 



2e, 



e+ X e_ 



(63) 
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Equations (|27l) and (1281) therefore reduce to 

= To^oj + (I)COJ' + J2 (e-b)C''-'", (64) 

b—x,y,z 



+ ((e-a)(e.b))C''^ (65) 

b—x,y,z 

or, explicitly, 

[1 - (1)] C"^' - (e,)C"J' - {ey)Cy' = ro(5oj, (66) 




FIG. 9: (Color online) Schematic view of a 3D topological 
insulator with a surface state subject to (a) perpendicular and 
(b) parallel magnetic field B. In the latter case the B- field 
dependence of the surface weak-antilocalization conductivity 
is due to the magnetic flux through the effective thickness of 
the surface state, determined by the decay length, A, into the 
bulk. 



+ [1 - (e^)] C-^- = ro<5,„ (67) 



{ey)C'^ + [l-{el)]Cy= =ro5y„ (68) 



[l-(ei>]C^^=ro<5,,. (69) 

In the above equations we use the short-hand notation 
e = e+ for the unit vector ekpn whose in- and out- 
of-plane components are given by Eq. (|21l) . Solving 
Eqs. - for the diagonal Cooperon coefficients, 
we have 



C""(q,c.) = 



1 



Dq2 



(Jxx 



Cvy 



D = 

„2 



(70) 



To 



1 



(ei>- 



^*n,wC0"(q,c.),(71) 
(72) 



Note that the singlet Cooperon C*'" ([70)1 remains gap- 
less also in the presence of the warping because it does 
not break the time-reversal symmetry. The warping only 
modifies the diffusion constant D through the transport 
scattering time (|62p . The triplet Cooperons C^^ and C^^ 
are already averaged, for presentation purposes, over the 
directions of q. In the absence of the symmetry breaking 
[cf. Eqs. dH]) and (g!])] and for Ttr/r^ < 1, the triplet 
Cooperons C^^ and C^^ as well as C^^ can be neglected 
compared to C*^'' in the conductivity corrections. 



D. Hikami boxes and net conductivity correction 



e^Dr 



dq 



[C°° - C^^] ,(74) 



Keeping only the singlet Cooperon C°°, we obtain the 
following expression for the net correction: 



e^D f dq 



XX 

2 



(|^C™(q)(75) 
(76) 



After the partial cancellation in Eq. (1751) : Dr{Ttr/T)'^{l — 
e^/2) = D {Ttr/r) {r/rtr) — D, the prefactor in the final 
equation ([7S)) depends only on the transport scattering 
time through the diffusion constant D [see, Eq. ((70)) and 
Ref. [4l| . Inserting Eq. (1701) into Eq. (|76l) , integrating over 
q with the upper cut-off (Dt)^^/^, and replacing — iw — >■ 
T^^, we obtain Eq. (|4]), already discussed in Sec. ID 

Below we will focus on the dependence of Saxx on the 
strength and oriention of an external magnetic field B, 
which can be used in experiments to identify the surface 
states in three-dimensional topological insulators. 



E. Magnetoconductivity in perpendicular field 

For B applied perpendicularly to the surface (see Fig. 
|9^) we write Eq. (ITO)) for the singlet Cooperon in the 
two-dimensional position representation^ 



Repeating the calculations described in subsec- 
tion nil Dl we express the conductivity corrections Sai]] 
(05]) and Saxx (SSI) in terms of the diagonal Cooperons: 



2(nlel)C- 



dq 



[2( 



)C 



00 



(2^)2 

2{nlel)Cyy-2{nlel)C^ 



1,(73) 



D[i\/ + yA(r) 



C°°(r,r') -(5(r-r'),(77) 



including the vector potential A(r) = (—By, 0,0) of the 
magnetic field (r = (a;, y, 0)). The solution is given by the 
Hilbert-Schmidt expansion in the Landau wave-function 
basis, which yields the well-known expression^^ for the 
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magnetoconductivity Acrxx{B) = 5axx{B) — 5axx{^)'- 



A<Jxx{B) 



e 

2Trh 



B B 



2el 







B 


J- surface 




B 


surface 



where the magnetic field B± corresponds to the 
Aharonov-Bohm flux of order oi h/e through a typi- 
cal area encircled by the interfering trajectories,— itr = 



Vp Ttr and 



T^, and Tpix) is the digamma function. 



FIG. 10: (Color online) Quantum-interference magnetocon- 
ductivity /S.axx{B) [in units of /2-Kh, see Eqs. (|78[) and 



for in- and out-of-plane field orientations. The mag- 
netic field scales, on which A(T2,i(_B) decreases, are fixed such 
that B /Bi_ = 2^/1^/ \ = 20. 



F. Magnetoconductivity in parallel field 

In the case of the parallel magnetic field B (see Fig. [Hb) 
the vector potential, A(z) = (B x z)z, depends explicitly 
on the coordinate z perpendicular to the surface (z is the 
unit vector). Therefore, the penetration of the surface 
state into the bulk needs to be taken into account. To 
do so we first transform the diffusion operator in equa- 
tion [Dc^ — icL>]C'^"(q) = 1 into the three-dimensional 
position representation and then make the Peierls sub- 
stitution iVr zVr -I- ^A(z) as follows 

/ drdze-''^'^x*{z)[D{iVr + f ^(2))^ - iw]e*i''x(^) 
xCOO(q) = l, (80) 

where the in-plane (r) integration goes over the surface 
area a, and the out-of-plane (2) integral involves the nor- 
malized wave function of the surface state, x(-2^)j which 
decays exponentially for z — >■ — c» on the length A in- 
versely proportional to the bulk band gap)2^ 

In the latter form (i.e. written with |z|) this equation can 
formally be used in the entire space — c» < z < 00. This 
helps to simplify further calculations because the system 
is symmetrically extended to the other half-space, < 
2; < 00, and the z integral in Eq. ([50]) can be calculated as 
J dz... = (1/2) J^^dz.... Having done this integration, 
we return to q representation of the Cooperon: 

C°"(q,c.)^ ^ . , (82) 

Dq^ + Tg — luj 

= 2De^B^\^/h^ = DX'^/2e%, (83) 

where tb is the time-scale for the suppression of the quan- 
tum interference by the magnetic flux through the thick- 
ness of the surface state, A. The latter is assumed much 
smaller than the magnetic length £b — \/^l'2\eB\. Note 
that Eq. (j83|) has a six-time larger numerical coefficcnt 
than the result for the quasi-2D quantum wellsj^ which 
reflects the difference in the electron confinement. 



Using Eqs. (f76|) and (f82|) we calculate the magneto- 
conductivity t^OxxiB) = Saxx{B) - SaxxiO): 

A...(i.).-^ln(l + ^) 

B, = 2^^B^. (85) 

Clearly, for a sufficiently small penetration length A <C 
2^ Itriip the in- and out-of-plane geometries have dis- 
tinctly different magnetic- field scales, B^^ ^ on 
which t^Oxx (B) decreases with B [see Fig. [TU] . The in- 
plane magnetoconductivity Aaxx{B) ([84]) can be veri- 
fied against recent magnetotransport measurements on 
BiaTea (Ref. HI. 
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Appendix A: Derivation of Cooperon equations (|25p 

The equation for the coefficients C'^ = 
E„s<..«. follows from Eq. 

X ^::pG2y (k, e + hw)G;,, (q - k, e)%,s] C'^' (q) • 
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It can be rewritten in a simpler form with the help of the 
identity 

J2 "^"pG^y (k, e + hiu)G%, (q - k, e)%,s, = 



Tr [g"""^ (q - k, e)*'1"G" (k, e + ntj)^'' 



(A2) 



where Tr and T denote the trace and transposition oper- 
ations, respectively, in a space. We therefore have 



ttNt, 







xG (q-k,e)r'^G (k,e + ;iw)«''' G'^^(q).(A3) 
For the orthonormal basis functions given by 



(J 



(A4) 



Eq. (jA3p assumes the following form: 

-2 n 



To - 2'nNTn 
G-* (q - k, e)cr*G" (k, e + huj)a' \ C'^ (q). (A5) 

J 



Here the tilde denotes the operation 



G^ = cr^ G 



(A6) 



which flips the pseudospin: or = cr^er^cr'^ = — <t. Thus, 
for q Eq. (jA5l) describes interference between the 
state with k, cr and its " time-reversed" partner with 
— k, — cr. 



Next, we evaluate the k integral in Eq. (IA5[) using the 
sharpness of the Green's functions at the Fermi level un- 
der condition (fT4|) and the smallness of the Cooperon 
momentum q and frequency w [see, Eq. ([35])]. In the de- 
nominator of G'*q_k it is sufficient to keep only the linear 
term in q, i.e. ^fq-k ~ ^ ^vq + ... (which should be 
compared with energy fi./T): 



q-k 



1 cr" — (Te_k 

2 - + fi^vq' 



(A7) 



At the same time, in the numerator we ap- 
proximate eq_k ~ e_k, neglecting small terms 

AqUlA^kl +Ml^ ^ Aq/E, « 1. 



With Eqs. ^9^, (TAT)) and under condition (|T4]) the k 
integral can be evaluated as follows 



^G-(q - k, .).^G^(k, . n.).^ -U'-P^f ^(Ck, n).a ^f'7-:y^";+'^^,-^f (A8) 



T / ^ cr -(Te_fc^.n V cr' -t-(Tefc^.„ a / da 77 ^ 777 7 

4 7 27r 7 (a - - ftw^n • q)(Ck - e„) 

4 J 27r z 



h/r + huj — hvpii ■ q 



27rA^T 1 



2tt 



(cr" - cre_fc^.„)(7*(cr" +o-efc^.„)cr^ 
1 — irw + irVpii ■ q 



(A9) 
(AlO) 

(All) 



where the DOS A^(a, n) is replaced by its Fermi surface value N and, then, the ^ integral is calculated in the complex 
plane. Inserting Eq. (jAlip into Eq. (K5\i yields Eq. ((25| . 



Appendix B: Evaluation of k integrals in Hikami 
boxes in Eqs. (|43]) and ([il]) 

To calculate the k integral in Eq. we first expand 
the Green's functions G^_j^ ^_^_f^^ and G^_j^ ^ in small 
Cooperon momentum q as we did in Eq. (jA7| of Ap- 
pendix |X| 



G 



cr" + cre_k 



(Bl) 



G, 



q— k,€ 



1 0-0 + cre_k 

2 - a + ^vq 



(B2) 



Then, the integral is converted to J J Nd^k--- and 
the S^k integration is again done in the complex plane 
following the same steps as in Eqs. (lASp - (jAlip . The 
final result is 
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where (...) = /g'^ ...(i(/)n/27r is the averagmg over the mo- 
mentum direction n. As the expression above is indepen- 
dent of e, the energy integral in Eq. p3)) is fi,/27r, which 
along with Eq. dBH]) yields Eq. (gS]). 



To evaluate the integrals over k and k' in Eq. pij) we 
set q, — in the single-particle Green's functions and 
re-group them as follows 



J (27r)2 J (27r)2 (^k,e "^^k Gk,c+;ii.jGq-k',e+rii.j)7'^(^q-k,c+;iw ^k'^c+Rcj "l^k' Gk',c)7^' ~ (^4) 

« / (^[G'(k)V-(k)G"(k)]y^,G;,^(-k) 1 ^G;^^(-k')[G"(k')V-(k')G^k')]7./5'. 
Each of the integrals can now be done in the similar way as in Eqs. (jASp - (IA11[) of Appendix \X\ 



/ (^[G'(k)V-(k)G"(k)],,^,G;,^(-k) 



ft2 23 



X [(aO + (Tefe,.„)V^^.„K-f *Tefe,.„)]y^,K + ae_fe,.„)^^^, (B5) 
G;,,(-k')[G^(k')V^(k')G^(k')]7..^ ' 



Inserting these into Eq. pij) we obtain Eq. 
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